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NONLINEAR  LUMPED  PARAMETER  MATHEMATICAL  MODEL 
OF  DYNAMIC  RESPONSE  OF  THE  HUMAN  BODY 

Gordon  R.  Hopkins 

West  Virginia  University 


ABSTRACT 

Two  nonlinear  models  of  man's  dynamic  response  to  lov  frequency  vi¬ 
bration  are  discussed.  The  first  model  uses  linear  spring  and  damper 
elements  but  accounts  for  the  nonlinear  geometry  of  visceral  mass  motion. 
This  model  adequately  reproduces  both  the  input  mechanical  impedance  and 
vibration  transmission  characteristics  for  a  seated  human  subject. 

The  second  model  includes  the  nonlinear  effects  of  the  lungs.  The 
influence  of  this  nonlinearity  on  the  dynamic  response  is  discussed  and 
compared  to  experimental  results  from  tests  on  animals. 

INTRODUCTION 

Over  the  past  two  decades  considerable  attention  has  been  paid  to 
the  effects  of  vibration  and  impact  on  the  human  body.  As  humans  are 
exposed  to  the  implements  01’  modern  technology  these  mechanical  stimuli 
reach  magnitudes  that  can  cause  physiological  damage  ranging  from  mild 
discomfort  to  death. 

The  inaugural  research  effort  in  the  study  of  man's  dynamic  response 
was  to  establish  human  tolerance  limits  to  vibration  and  impact.  Once 
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these  vere  established,  the  attention  of  investigators  turned  toward  the 
bore  perplexing  problem  of  understanding  the  injury  mechanisms  that  set 
the  tolerance. 

Using  human  subjects  as  test  specimens  to  gather  information  about 
the  internal  dynamics  of  the  body  has  not  always  been  a  gratifying  expe¬ 
rience,  Yet,  animal  substitutes  in  most  cases  have  proved  equally  dis¬ 
appointing.  The  problems  associated  with  trying  to  glean  information 
about  the  dynamics  of  the  internal  organs  and  structural  components  of 
the  human  body  is  complicated  by  the  necessity  of  having  to  gather  data 
without  damage  to  the  subject.  This  limits  the  types  of  instruments  that 
may  be  used  as  well  as  severity  of  the  test  environment.  In  effect,  in¬ 
vestigators  have  been  gathering  data  about  internal  dynamics  of  the  human 
body  by  sensing  only  their  external  effects. 

Useful  data  have  been  gathered  in  this  manner  however.  For  instance, 
most  of  our  understanding  of  the  body's  dynamic  response  to  low  frequency 
vibration  has  been  based  on  measurements  such  as  driving  point  mechanical 
impedance ,  external  strain,  vibration  transmission,  and  cineradiographs. 

In  order  to  make  the  best  use  cf  these  data,  mathematical  models  have 
often  been  employed  in  the  effort  to  expand  our  understanding  of  the 
human  body  as  a  dynamical  system. 

This  report  discusses  one  such  model.  A  nonlinear  lumped  parameter 
model  that  explains  phenomena  observed  during  low  frequency  vibration  of 
a  seated  man,  such  as  the  second  and  third  spinal  resonances. 
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BACKGROUND  INFORMATION 

The  rationale  of  the  model  is  best  discussed  by  first  briefly  re¬ 
viewing  the  experimental  data  on  which  it  is  based. 

Mechanical  impedance  has  proved  to  be  a  valuable  and  well  used  tool 
in  studying  the  vibration  response  of  biological  systems.  With  its  use 
as  a  parameter,  the  dynamic  response  of  a  human  subject  to  vibration  can 
be  studied  as  a  "black  box"  much  as  one  would  an  unknown  electrical  cir¬ 
cuit  with  input  electrical  impedance.  Although  it  has  limitations,  this 
approach  yields  much  information  about  man's  dynamic  response  to  sinus¬ 
oidal  vibrations  by  showing  system  resonances. 
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Fig.  1  Mechanical  Impedance  and  Phase  Angle  for  Seated 
Subject ^ ^ 
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illustrates  the  mech- 


1  2 

Figure  1  based  on  the  work  of  Coermann  et  al,  * 
anical  impedance  of  a  typical  subject  to  lew  frequency  sinusoidal  vi¬ 
bration.  The  predominant  features  of  this  curve  are  the  peaks  at  5*  11, 
and  lU  Hz.  These  peaks  suggest  points  of  resonance.  As  von  Gierke  sug¬ 
gests,  the  first  peak  is  due  to  the  abdominal  resonance,  and  the  second 
and  third  due  to  some  action  of  the  upper  torso  on  the  vertebral  column. 

A  second,  and  equally  useful  tool  used  in  vibration  analysis  is  vi¬ 
bration  transmission,  or  transmissibility.  Again  Coermann1  shows  us,  and 
is  substantiated  by  the  work  of  Pradko  who  used  random  vibration  instead 
of  sinusoidal,  that  the  vibration  response  of  a  seated  human  subject  can 
be  characterized  by  transmissibility.  Figure  2  shows  this  result. 


0  5  10  15 

Frequency,  Hz. 


Fig.  2  The  Transmission  of  Vibrations  from  the  Seat  to  the 
Head  of  a  Seated  Human  Subject'^-' 
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The  transmission  factor  for  a  seated  subject  shows  three  peaks  nominally 
at  5*  11,  and  l4  Hz.  As  vith  the  mechanical  impedance,  it  is  hypothe¬ 
sized  that  these  peaks  are  the  result  of  resonant  systems  in  the  body. 
One  should  note  that  the  peak  corresponding  to  the  abdominal  resonance 
transmits  the  largest  amplitude  vibrations  to  the  head. 

Another  noteworthy  set  of  parameters  that  characterize  the  body's 
dynamic  response  to  low  frequency  vibration  are  the  strain  measurements 
made  by  Clark  et  al.^  and  shown  in  Figure  3.  It  is  observed  from  these 
measurements  that  the  first  peak  for  the  upper  and  lower  abdominal  strain 
occur  at  approximately  1  Hz  apart.  The  implication  here  is  that  the  vis¬ 
ceral  mass  works  against  a  different  spring  when  it  moves  toward  the  head 
than  it  does  when  it  moves  toward  the  pelvis. 


Frequency,  Hz, 


Fig.  3  Strain  for  Seated  Subject  During  Low  Frequency 

Vibration^) 
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The  final  set  of  measurements  to  be  considered  are  those  made  by 
White  et  al,^  of  colon  pressure  during  low  frequency  vibration.  As  anti¬ 
cipated,  these  measurements  show  the  effect  of  the  abdominal  resonance  in 
agreement  vith  the  lower  abdominal  strain  shown  in  Figure  3.  If  the  sec¬ 
ond  and  third  peaks  of  Figures  1  and  2  are  due  to  vertebral  column  reson¬ 
ances  as  hypothesized,  then  the  effect  on  colon  pressure  would  be  minimal. 
As  can  be  seen  in  Figure  4,  this  is  precisely  the  case. 


Fig.  L  The  “Mean"  Colon  Pressure  During. Low  Frequency 
Vibration  for  a  Seated  Subject'®' 

In  summary,  th«se  experimental  measurements  require  a  model  toi 

1.  Have  a  similar  impedance  and  phase  versus  frequency  re¬ 
lationship  as  that  for  a  human  subject, 

2.  Demonstrate  the  same  vibration  transmission  characteristics 
as  does  a  human  subject* 
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3,  Corroborate  the  organ  motions  as  inferred  from  the 
measurements  of  strain  and  colon  pressure* 

DEFINITION  OF  MODEL 

It  has  been  demonstrated  by  Hopkins J  Suggsf  and  others  by  the  con¬ 
struction  of  mathematical,  electric  analog,  and  mechanical  models  that  an 
uncoupled  two  degree  of  freedom  system  as  shown  schematically  in  Figure  5 
can  adequately  simulate  the  input  mechanical  impedance  of  the  human  body 
over  the  0-15  Hz  range.  For  purposes  where  only  the  input  impedance  is 
necessary,  such  as  coupling  with  other  machinery  for  their  evaluation, 
this  type  of  model  is  sufficient.  However,  for  inquiring  into  the  inter¬ 
nal  dynamics  of  the  human  body,  such  a  model  is  not  satisfactory.  It 
lacks  sufficient  anatomical  similarity  to  predict  injury. 


Fig.  5  Uncoupled  Two  Degree  of  Freedom  Model 
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There  is  no  vay  to  use  this  type  of  model  to  study  the  transmission 
of  vibration  to  various  parts  of  the  body  or  to  corroborate  such  data  as 
transmissibility,  strain  and  internal  pressures.  This  deficiency  leads 
one  to  consider  a  coupled  mass  linear  model^ such  as  shovn  in  Figure  6 , 
that  could  accommodate  these  considerations. 


Fig.  6  Linear  Coupled  Mass  Model 

The  linear  coupled  mass  offers  enough  geometric  similarity  to  the  human 
body  to  be  able  to  relate  components  to  their  anatomical  counterparts. 

For  instance*  and  could  represent  the  vertebral  column,  Mg  the 
upper  torso  mass,  the  visceral  mass,  the  lover  torso  mass  and  so  on. 
However,  this  model  does  not  adequately  simulate  the  input  impedance. 
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Vertebral 

‘Column 


Fig.  ?  nonlinear  Geometry  Model  of  Uun&n  Body 


If  one  examines  the  anatomy  of  the  human  body  in  more  detail,  it  may 


be  observed  that  the  visceral  organs  are  not  tethered  but  are  supported 
in  the  abdominal  cavity  by  the  abdominal  muscle  wall,  the  pelvis,  and  the 
diaphragm.  When  the  visceral  mass  moves  during  vibration  it  exhibits  a 
nonlinear  geometry.  As  the  visceral  mass  moves  upward  toward  the  head, 
it  compresses  the  lungs  but  does  not  put  tension  on  the  abdominal  wall  or 
pelvis  since  it  is  not  tethered  to  either.  On  the  other  hand,  as  the  mass 
moves  toward  the  pelvis ,  it  does  not  put  tension  on  the  diaphragm  and  in 
turn  the  lungs. 

The  simplest  model  that  can  accomodate  this  nonlinear  geometry  ia 
one  with  linear  elements  constructed  such  that  the  springs  are  not  rigid- 
ly  attached  to  the  mass  that  represents  the  visceral  mass. 

This  model  is  shown  schematically  in  Figure  7  where  the  linear 
spring,  k^,  and  damper,  C^,  represent  the  vertebral  column,  kg,  the  ab¬ 
dominal  wall,  k^,  the  lungs  and  Cg,  the  frictional  damping  on  the  inter¬ 
nal  movement  of  the  organs, 

DWAHIC  BEHAVIOR  OF  NGSLX8EAR  GEOMETRIC  MODEL 

A  mathematical  model  based  on  the  schematic  shown  in  Figure  7  has 
been  fo resale ted  and  interpreted  as  a  computer  model  using  IBM's  Contin¬ 
uous  System  Modeling  Prcgras.  The  dynamic  bebavf  or  of  this  model  is 
shown  in  Figures  8,  9,  and  10. 
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DISCUSSION  OF  NONLINEAR  GEOMETRY  MODEL 
Comparison  of  Figures  1  and  8  shows  that  this  two  degree  of  freedom 
model  with  nonlinear  geometry  adequately  reproduces  the  mechanical  im¬ 
pedance  and  phase  angle  characteristics  of  the  human  body}  even  to  dis¬ 
playing  the  second  and  third  resonant  peaks.  Comparison  of  Figures  2  and 
9  shows  that  the  model  also  demonstrates  the  same  general  vibration 
transmission  characteristics  as  the  human  body. 

If  one  assumes  that  the  displacements  of  the  visceral  mass  can  be 
related  '■>  the  external  strains  on  the  chest  and  abdominal  wall,  then 
Figure  10  agrees  with  the  strain  measurements  presented  in  Figure  3, 

It  is  concluded  that  the  model  does  satisfactorily  demonstrate  the 
same  dynamic  response  to  low  frequency  vibration  as  does  a  seated  human 
for  small  amplitude  vibrations. 

NONLINEAR  MODEL 

Since  all  the  data  on  which  the  above  model  is  based  was  gathered 

from  forced  sinusoidal  vibration  with  low  amplitude  accelerations,  from 

lA  to  1/2  g. ,  Its  use  for  the  study  of  vibration  with  higher  amplitude 

acceleration  or  impact  is  questionable.  Experience  with  the  nonlinear 

geometry  model  has  shown  that  its  dynamic  response  is  not  dependent  on 

the  amplitude  of  the  forced  vibration, 

o 

Krause  and  Lange  have  shown  with  experiments  on  pigs  that  mech¬ 
anical  impedance  is  nonlinear  with  respect  to  the  acceleration  amplitude 
of  the  forced  vibration.  Their  results  are  reproduced  in  Figure  12, 
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Fig.  12  Impedance  of  70- lb  Pig^ 

If  the  nonlinear  geometry  model  is  slightly  sophisticated  tc  include 
the  nonlinear  effects  of  the  lungs  by  modeling  them  as  a  piston  in  a  cyl¬ 
inder  with  an  orifice,  then  the  nonlinear  behavior  of  biological  systems 
observed  by  Krause  can  be  simulated. 

The  equations  of  motion  for  this  model^  shown  in  Figure  11,  are: 

■  hjUg  -  +  CjUg  -  x^  +  kg(x3  -  x^)  ♦  C3(x3  -  x^) 

♦  6(*3  -  x^  +  F(t) 

®2*2  *  -k(xg  -  xL)  -  C1(xg  -  xx)  + 
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»2*3  -  **3  « ^  -  xx)  -  C3(x3  -  »j.)  -  F] 


LUNGS 


where 


«(x) 


V/ 


x  >  0 


X  B  0 


The  force  exerted  by  the  lunge  on  the  visceral  mass  (F^gg)  is  de¬ 
rived  from  the  model  shown  in  Figure  1^« 


Total  Volume 
(366  in3) 


Fig.  A**  Lung  Model 

Assuming  isotropic  compression,  the  equation  of  state  pV  »  mRT  gives  the 
following  formulation  for  F^^ggS 
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P 


J[c[2]8[i|  -U-l)|)-4} 


dt 


where 


V  ■  A^fL  -  (x3  -  x2)]  , 

-\  f  nf  Pa  2/ 6  p  d— 1 

-  -  AtVrrr  pL  >l  «<f 1  [1  -  <r>'r  J 


where  s  denotes  smaller  of  p  and  p 


and  L  the  larger  of  the  two, 
and 

r. 

dt 


*  •  "Wtii  -  J " 

Then, 


F  s 

'lungs 


tU2  -  *3>  , 

(Xj  -  X  )  p  \ 


sign  (p  -  pL) 


With  this  mathematical  formulation  a  digital  computer  model  was  con¬ 
structed  using  IBM’s  Continuous  System  Modeling  Program. 

The  CSMP  for  this  nonlinear  model  was  subsequently  used  to  deter¬ 
mine  the  effects  of  different  amplitude  levels  of  vibration  on  the  mech¬ 
anical  impedance.  Figure  15  shows  these  results. 

Although  there  are  no  mechanical  impedance  data  available  for  hu¬ 
mans  at  low  frequency  vibration  levels  as  high  as  3  g's.,  comparison 
with  the  results  of  Krause,  Figure  12,  shows  that  the  dynamic  behavior  of 
the  nonlinear  model  is  similar  to  that  of  a  pig.  That  is,  there  is  a 
general  stiffening  of  the  system  as  the  acceleration  level  is  increased. 
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CONCLUSIONS 


1,  The  inclusion  of  nonlinear  geometry  in  the  two  degree  of 
freedom  model  of  the  human  "body's  dvnamic  response  of  the 
human  body  improves  the  dynamic  similarity, 

2,  The  nonlinear  mechanics  of  the  lungs  probably  account  for 
most  of  the  body's  nonlinear  response  to  low  frequency 
vibration. 

3,  If  sufficient  effort  were  expended  to  establish  the  gov¬ 
erning  parameters  of  the  nonlinearity  by  vibration  tests, 
nonlinear  lunged  parameter  models  could  add  understanding 
of  the  body  mechanics  during  certain  impact  situations. 
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